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Abstract 

A series of recent articles introduced a method to construct stochastic partial differential 
equations (SPDEs) which are invariant with respect to the distribution of a given condi- 
tioned diffusion. These works are restricted to the case of elliptic diffusions where the drift 
has a gradient structure, and the resulting SPDE is of second order parabolic type. 

The present article extends this methodology to allow the construction of SPDEs which 
are invariant with respect to the distribution of a class of hypoelliptic diffusion processes, 
subject to a bridge conditioning, leading to SPDEs which are of fourth order parabolic type. 
This allows the treatment of more realistic physical models, for example one can use the 
resulting SPDE to study transitions between meta-stable states in mechanical systems with 
friction and noise. In this situation the restriction of the drift being a gradient can also be 
lifted. 



1 Introduction 

In previous works, see e.g. |SVW04[ IHSVW051 IHSV07) or |HSV09] for a review, we described 
an SPDE-based method to sample paths from SDEs of the form 

x{t)^ f{x{t)) +w{t) Vie[0,T] (1) 

where w is white noise, conditioned on several different types of events. The method works by 
introducing an 'algorithmic time' r and constructing a second order SPDE of the form 

drX{T,t)=d^x{T,t)+J\f{x{T,t))+V2drW{T,t) V(t, t) G R+ X [0, T] (2) 

which has t as its space variable. Here drw{T,t) is space-time white noise. The nonlincarity 
Af and the boundary conditions of the differential operator are constructed such that, in 
stationarity, the distribution of the random function t i— >■ x{t, t) coincides with the required 
conditioned distribution. (See also |RVE05 '.') It transpires that the distribution of ^ under the 
bridge conditions x(0) = x{T) = corresponds the choice 

M,{x) = -U{x)d,h{x) -\dlMx) 

(written using Einstein's summation convention) and use of Dirichlct boundary conditions for &f . 

Assuming ergodicity of the sampling SPDE, one can now solve the sampling problem by 
simulating a solution to Q up to a large time r and then taking <: i— >■ x(t, t) as an approximation 
to a path from the conditioned SDE. The resulting sampling method has many applications, 
some of which are described in [AHSVOT] . The biggest restrictions of this method are that the 
derivation requires the drift / to have some gradient structure and the diffusion matrix (chosen 
to be the identity matrix in ([Ij above) to be invertible. 

In this article we consider the different problem of sampling conditioned paths of the second 
order SDE 

mi{t) ^ f {x{t)) ~ x{t) + w{t) vte[o,r], (3) 
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conditioned on x{0) — X- and a;(T) — xj^. Equation ^ could, for example, describe the time 
evolution of a noisy mechanical system with inertia and friction. Rewriting this second order 
SDE as a system of first order SDKs for x and x leads to a drift which is in general not a gradient 
(even in the case when / itself is one) and, since the noise only acts on i, one obtains a singular 
diffusion matrix. Thus, this problem is outside the scope of the previous results. However, it has 
enough structure so that it still can be treated within a similar framework. Indeed, we derive in 
this article a fourth order SPDE of the form 

d^x(T,t) = {pi -n?df)x{T,t) + Af{x){T,t) + V2drw{T,t) V(T,t) e M+ X [0,T] (4) 

where, again, the boundary conditions and the drift term M are chosen in such a way that the 
conditioned distribution of ^ is stationary for (U). 

One surprising fact about this result is that it does not require / to be a gradient. In our 
earlier works, even the appropriate definition of solutions for the (formal) second order SPDE 
derived to sample conditioned paths of ^ in the non-gradient case is not clear (see |HSV07[ 
section 9] or |AHSV07l section 9.2]); the analysis for elliptic equations is thus restricted to 
the gradient case. In contrast, the greater regularity of solutions to SPDE here, sampling 
conditioned paths of ([3]), allows us to obtain existence results for the fourth order SPDEs arising, 
without any gradient requirements on /. 

In the special case where / is a gradient and /(x„) = fi^+) = Oj the components of the 
nonlinearity J\f can be written as 

Afjix) = ~fi{x)djfi{x) -mdtXt dtXk dj^f.^x) + mdt(dtXi{difj{x) + djf.i{x))^ 

(using Einstein's summation convention again). It is tempting to try to derive ^ by taking the 
limit m J, in ([U , in particular since the first terms of the corresponding nonlinearities coincide. 
It transpires that taking this limit is not entirely trivial: one needs to argue that on one hand, 
mdtXidtXk — >■ \5ik as m J, 0, but that the term mdt{dtXi{difj + djfi)) becomes negligible in 
the limit. Nevertheless, this argument can be made exact; see |HailOj for a rigorous derivation 
of the required limiting procedure. 

One novelty of this article compared to earlier work like jZab88[ IHSV07| is that there is no 
natural Banach space (like the space of continuous functions) on which the nonlinearity is well- 
defined and on which the linearised equation generates a contraction semigroup. The reason for 
this is that the linear operator of the equations studied in this article is a fourth-order differential 
operator. Another technical difficulty stems from the fact that the nonlinearity Af has very weak 
dissipativity and regularity properties. 

While preparing this text, we performed some numerical simulations on the fourth order 
SPDE presented here. Our aim was to study whether the SPDE could be used as the basis 
of an infinite dimensional MCMC method. Differently from the situation in earlier articles, 
these simulations proved prohibitively slow and the resulting method does not seem like a useful 
approach to sampling. This is mainly due to the fact the convergence time to equilibrium seems 
to grow like and thus can get very big for non-trivial problems. In the gradient case, since the 
system converges to the second-order SPDE as m — 0, one could expect improved convergence 
rates for small values of m. However, the theory developed in [HailO] suggests that the relevant 
lengthscale for the small-m problem is m, suggesting that one would need numerical simulations 
that resolve significantly smaller scales than that in order to obtain reliable results. Again, this 
would lead to inefficient numerical methods even in the case of small m. Consequently, we do 
not include our simulation results in this article. 

For a number of articles considering fourth-order (S)PDEs, see for example [BMPWOTj . 
[DPD96| and |LCM96| . Alternative methods to construct solutions of SPDEs and to identify 
their stationary distributions are based on the theory of Dirichlet forms, see e.g. [RM92) . 

The text is structured as follows: in section [5] we give a detailed description of the sampling 
problem under consideration and formulate the main result in theorem 01 The proof of this 
result is given in sections |3l |4] and [5] 

Notation 

Throughout the article, we will use the notation as introduced above: by G [0,r] we denote 
'physical time', i.e. the time variable in equations like ([ij and ^ which define the target 
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distributions. By ct, r > we denote 'algorithmic time', i.e. the time variable in sampling 
equations like ^ and (|3]). Thus, in the sampling SPDEs, r takes the role of time and t takes 
the role of space. 
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2 The Sampling Problem 

In this section we give the full statement of the sampling problem we want to solve; the main 
result is contained in theorem 21 

First consider the following unconditioned second order SDE: 

mx{t)^ f{x{t))^x{t)+w{t) Vie[o,r], 
a;(0) = xo, i(0) = wq 

where the solution x takes values in R'', m > is a constant, / : M'' M'' is a given function, and 
w is a standard Brownian motion on . The initial conditions xq and wq are either deterministic, 
or random variables independent of w. The solution to this SDE can be interpreted as the time 
evolution of the state of a mechanical system with friction under the influence of noise. In 
this case m would be the mass and / would be an external force field. Models like this are, 
for example, widely used in molecular dynamics since for conservative forces / they describe 
Hamiltonian systems in contact with a heat bath. In this context, equation ([5|) is called the 
Langevin equation. The limiting case m = corresponds to the Brownian dynamics ([T|). 

Remark 1. Arbitrary constants in front of the x and w terms can be introduced using a scaling 
argument: Let /3,7 > and define the process y by y(t) — x(t/^)/ •\//3/2. Then y solves the 
SDE 

mm^f{y(t))-"fyit) + \lj'^ii) Vie [o,f] 

where rh = 7^771, f{x) — /(\//5/2a;) / y/ (3/2 and T = 7T. Thus, by rescaling T, m and F we can 
assume (3 = 2 and 7=1 without loss of generality. 

For our analysis we rewrite the second order SDE ([5]) as a system of first-order SDEs in the 
variables x and x. We get 

dx{t) — x{t)dt, x{0) = xq, 

mdx{t) ^ f{x{t))dt~x{t)dt + dw{t) x{0) = vq. 

In the Hamiltonian case f{x) — ~WV{x), and provided that the potential V is sufficiently 
regular, it can be checked that the Boltzmann-Gibbs distribution 

exp^— 2(y(a:;) + -mx^)^ dxdx 

is invariant for ([5]) . If V is sufficiently coercive, this distribution can be normalised to a probab- 
ility distribution. Note that in equilibrium, the position x and the velocity x are independent. 
Thus, in stationarity, the velocity satisfies x{t) ^ Af{0, l/2m) for all t e [0, T]. We will, even for 
the non-gradient case, use this distribution for the initial condition for x. 

Definition 2. For T > 0, x^,x+ e and /: R'' -> R'', define Q"'^" to be the distribution of 
the process x given by © where a;o = X- and vq ^ A/'(0, l/2m), independent of w. Define the 
target distribution iq i^g ^i^^g distribution of x under Q^'^" , conditioned on x{T) = 
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The sampling problem considered in the rest of this article is to find a stochastic process with 
values in ([0, T] , R'') which has the target distribution •^•■^+ as its stationary distribution. 
Note that is just the distribution of a; and not of the pair {x, i) and thus is a probability 

measure on L^([0, T],R''). Considering this distribution is sufficient since for solutions of (O 
the initial condition a;(0) — a;_ allows to find a bijection between the paths x and x. If / is a 
gradient, the distribution ' coincides with the distribution of the process in stationarity, 

conditioned on x(0) = x_ and x{T) = x+. 

Definition 3. Let L denote the formal differential operator 

L = -m^df + 

and define £ to be this differential operator on the space ([0, T], K*^) , equipped with the 
domain 

V{C) = {x e H'' I x{0) ^ x{T) = 0, md'^xiO) = dtx{0), md^x{T) = -dtx(T)} (7) 

where ~ H'*([0,T],K'^) is the Sobolev space of functions with square integrable generalised 
derivatives up to the fourth order. Furthermore, let x : [0, T] -> R*^ be the solution of the 
boundary value problem Lx — with boundary conditions 

x{0) = X-, x(T) = x+, mdfx{0) = dtx{0), md^x{T) = -dtx(T). (8) 

We will see in lemma [T7] that the operator £ given by this definition is self-adjoint and 
negative definite. 

Theorem 4. Consider the ([0, T], M'*) -uaZitec? equation 

dx{T) = C{x{t) — x) dr +M{x{t)) dr + V2dw{T), x{0) = xq. (9) 

Here £ and x are given in definition]^ w is a cylindrical Wiener process, xq G L^([0, T], M'') 
and 

J^k{x) = ~fi{x)dkfi{x) + mdtXidtXj dfjfk{x) 

- dtXi{difk{x) - dkfiix)) + md^Xi(difk{x) + dkfi{x)) (10) 
+ m(fk{x-)dfSQ - fk{x+)dtdT) 

for k = 1, . . . ,d where we used Einstein's summation convention over repeated indices, Sq and 
St are the Dirac distributions at and T respectively, and all derivatives are taken in the 
distributional sense. Assume that f e C'^(]R'', K'^), that the partial derivatives dif and dijf are 
hounded and globally Lipschitz continuous for all i,j = l,...,d and that there are constants 
/? < 1 and c > such that \f{x)\ < \x\^ + c for all x G M''. Furthermore assume that the 
SDE ([6]) a.s. has a solution up to time T . Then the following statements hold: 

a) For every xq G L^([0,T],R'^), equation ^ has a unique, global, continuous mild solution 
with E(|j2:(T)||L2) < oo for all r > 0. 

b) For T > Q, the solution x{t) a.s. takes values in the Sobolev space H^([0,T],R'^). 

c) The distribution giyen hy definition]^ is invariant for 

Remark 5. The sub- linear growth-condition |/(a;)| < jxp -f c on the drift seems quite technical. 
The condition is only required for the bounds in lemmajH] We believe that an additional, linear 
drift term can be added by incorporating it into the linear operator £, following a similar 
procedure as in jHSVWOS"] . 

Remark 6. As in [HSV07[ remark 5.5], we see that the terms involving derivatives of Dirac 
distributions can be interpreted as modifications to the boundary conditions. Proceeding this 
way, we see that ^ is formally equivalent to the SPDE 

drXk(T,t) = Lxk(T,t) - fi{x)dkfi(x) +mdtXidtX.j dfjfk{x) 

- dtXi{difk{x) - dkfiix)) + md1xj{djfk{x) + dkfj{x)) + V2 drWk{T,t), 
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where drW is space-time white noise, endowed with the boundary conditions 

x{0) — X- = x^ 

mdfx{0) = dtxiO) + fix-), mdfxiT) = -dtx{T) + f{x+). 

In the one-dimensional case / : R M this SPDE simphfies fm'ther to 

drx{T, t) = Lx{t, t) - fix) fix) + m idtx)^ fix) + 2m d^x fix) + %/2 drwir, t). 

Remark 7. Using a standard bootstrapping argument hke the one in the proof of |Hai09[ 
Theorem 6.5] one can show that, in fact, the solution x of ([9]) takes values in the Sobolev space 
Wi[Q,T],W^) for every r < 3/2. 

The remainder of the article gives a proof of theorem 21 We start the argument, in section [31 
by collecting some results about the differential operator C. Section |31 shows that the theorem 
holds for the linear case / = 0, in which case A/" = 0. Finally, Section |5l then completes the 
proof by showing that introduction of the drift N changes the stationary distribution of ^ in 
the correct way to account for a non- vanishing /. 

3 Analysis of the Linear Operator 

This section collects some results about the operator C from definition [31 Since we are only 
interested in the operator itself and not in the full SPDE, in addition to the scaling-argument 
from remark [H we can rescale t. Thus, throughout section [3l we will consider the operator C 
defined as 

where 7 = , on the domain 

ViC) = |x e I x(0) = x(7r) = 0, — dUiO) = dtxiQ), — d^xin) = -dtxin)]. 

Then, after rescaling t, C differs from the operator C from definition |3l only by multiplication of 
a positive constant. 

Throughout the rest of the paper we will use the following notation: We denote by (S'(t))t->o 
the semigroup associated to C on H = L^([0, vr], R'') and by Ha = T^ii,^ £-)"') the associated 
interpolation spaces. 

3.1 Approximation to the Spectral Decomposition 

Lemma 8. C is a self-adjoint, negative definite operator on ([0, tt], R'') . 
Proof. Using partial integration it is easy to see that 

{x,Cy)~—vn? I dfxdfydt— / dtxdtydt 
Jo Jo 

- m(^dtxiO)dtyiO) + dtxin)dtyin)^ 

for all x,y € i.e. the operator C is symmetric and negative. Its self-adjointness can be 

checked just as in |RS72I section VIII] . ■ 

Lemma 9. Let Xk, fc e N &e the eigenvalues of —C and Cfc be the corresponding eigenf unctions. 
Define furthermore 

ifi''' I i = 1,2,3,4) = (sin/ci, coskt, e^*^*, e-''^"-'^). 

Then the following statements hold: 
a) The eigenvalues of C satisfy Xk = + Oik^). 
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b) There exist functions g'j^^ such that ek{t) ~ sin(fci) + ^ X]j=i 9k\t)fk^i^) '^^^ ^ ^ [Oj^ 



(i)/ 



and such that swpj^i sup^^gpj ||5["'''||c3 < for every j > 0. 

Proof. Since C acts independently on each coordinate, we can assume d 
generality. The eigenfunctions of C can be written in the form 



1 without loss of 



36" 



where 
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-^46" 



with \ = fj^ the corresponding eigenvalue. The coefficient vector ^ g is determined by the 
boundary conditions: for x to be an eigenfunction of £, ^ must satisfy ~ where 



1 

{aK+ - l)c-'^+'^ 
aK+ + 1 



1 



1 



1 



e — e 
aK+ + 1 — aK_ — i — aK_ + i 

(q!K+ — l)e^''+'^ {—aK^+i)e^^-'^ {~aK^ - i)e^'' 



\ 
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and a = T/ttj. Setting k — (k+ + k-)/2 = /i + ©(l/zi'^) for ease of notation, we note that this 
equation has non-zero solutions if and only if 

= detv4^ = 8i((Q!^K^ + an) sin(K_7r) - (| + an) cos(K_7r)) + 0(1//^) 
= 8za^/i^ sin fin + ©(/i). 

It follows immediately that, at least for large values of one has fi = k + 0{l/k) with fc S N, 
so that Afc = fc"* + ©(fc^) as requested. In particular, one has 



K± = k 



Pi 



0{\/k% 



for some constants j3± E M. It remains to check the statement about the eigenfunctions. 

Given that we already have good control on the eigenvalues, our claim will follow if we are 
able to show that one can choose ^ = (0, 0, 5,-5) + 0{l/k). Expanding Af^ in powers of k, we 
obtain 



A„ 



/O \ 



a —a —a 

—a —ctj 



/Oil 1 

1 (-l)*^ (-l)'^' 

1 -i i 

\1 ic{-l)'' -ic(-lfj 



+o(i/fc)^fc4°)+4i)+o(iA), 



for c — 1 ~ aj3^Tr. It now follows from standard perturbation theory (see e.g. [Kat66[ the- 
orem II. 5. 4]) that the eigenvector ^ with eigenvalue can be written as ^ = +0{l/k), where 
satisfies 4°^?^°^ = 0. Since is degenerate, this is however not sufficient to determine 
uniquely, but only tells us that is of the form {a + b,a + b, a, b) for a, 6 € JR. In order to 
determine a and 6, we have to consider the next order, which yields the compatibility condition 
A^Psy^^ G Range ^If''. This compatibility condition can be rewritten as a + 6 = 0, so that we 
can indeed choose — (0, 0, 5,-5), as requested. ■ 



3.2 The Relation Between Interpolation and Sobolev Spaces 

In this section, we show how the interpolation spaces Ha associated to the operator C relate to 
the usual fractional Sobolev spaces. These results are "well-known" in the folklore of the subject. 
However, in our specific context (especially since we need to consider fractional exponents), we 
were not able to derive them as straightforward corollaries from results in standard textbooks 
on function spaces, like |Tri83[ [Tri921 ITri06] . Because of this, and since one can find rather short 
and self-contained proofs, we prefer to include them here. 

Before we turn to this however, we start with a comparison between the interpolation spaces 
of the Dirichlet Laplacian and the periodic Laplacian. These are going to be useful in the sequel. 
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Let Aq denote the Laplacian on [0, tt] with Dirichlet boundary conditions and let A denote 
the Laplacian on [0, 27r] with periodic boundary conditions. These operators are self-adjoint in 
Hq = L^([0,7r]) and H — L^([0,27r]) respectively. We denote by Hg the domain of Aq^^ and 
by the domain of A*/^ (defined in the usual way through spectral decomposition). The aim 
of this section is to study the correspondence between these two different types of fractional 
Sobolev spaces. Denote by /, : Hg the map 



fit) for<G[0,^], 
-/(27r-t) for < e (7r,27r]. 



Note that t/2 is an isometry since it maps the eigenfunctions of Aq into eigenfunctions of A. 
This therefore defines an inclusion Hg C H'' . A natural left inverse for l is given by the restriction 
map 

^/ = /|[o,.] forall/eff. 

However, r is not an isometry and, for s > 1/2, it certainly doesn't map H'* into in general 
(since the constant function 1 belongs to every H" but only belongs to Hg for s < 1/2). We do 
however have the following: 

Lemma 10. The restriction operator r is bounded from H* into Hg for any s < 1/2. 

Proof. Note first that H* is isomorphic to H via the isomorphism x H> A'^^^x and similarly for 
Hg , so that the study of r as an operator from H* to Hg is equivalent to the study of the operator 
Ag^^rA"*/^ from H to Hq. Furthermore, we know that r is bounded from H to Hq, so that it 
suffices to show that the operator A ~ (Ag^'^rA^''/^ — r) is bounded from H to Hq. 

Since A maps sin(n • ) to for every n, it suffices to consider A on the subspace of H given 
by even functions, and generated by the basis of eigenfunctions of A given by (pnit) = ^ cos nt. 
Define furthermore the basis of eigenfunctions of Aq given by ipm = f sin mt. This yields for A 
the matrix elements 

2 

Aran = (V"™, Aip,^ = —{m^u^" - 1) / sm{mt) cos{nt) dt 



" JO 

2/ , — w T it m + rt IS odd, 

if m + n is even. 

It follows that there exists a constant C > such that ||^|| < C||A||, where the operator A is 
defined via its matrix elements by 

_ , if 71 > TO, 

n. 

Now it is a straightforward exercise in linear algebra to show that, given an orthonormal basis 
{<^n}n>Oj an operator A is bounded if there exist positive numbers /m.n such that the inequalities 

sup ^ \{A^Pn, Aipm)\ ^ ^ |(i(^„,l(/J„j)|/„.„j < OO, (11) 

">0^>g Jn,m ™>0„>o 

both hold. (Just expand for x — X]n>o^"'/^" ^^'^ make use of the inequality |a:;„Xm| < 

1_ ^m/ri.m ^jj]^ show that (|lll) does indeed hold for A as above. Assuming without 

loss of generality that m> n, we have the bound 

n m OO 

2s-2, 



fe=i 



< C(nTO"2 ^ m'-^n-') < Cm'-^ 



n 
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and similarly for n > m. Here, we have made use of the fact that s < ^ to ensure that the last 
sum converges. It remains to check that the bounds (|lip are satisfied for some choice of fm,n- 
With the choice frn,n = \/ fnjn^ we obtain 

ml 



oo 

3 1 



m>0 ' m=l 



where we made again use of the fact that s < ^. The second bound in is obtained in an 
identical way with the roles of m and n reversed. ■ 

For s > 1/2, the problem is that elements of Hg are forced to be equal to at the boundary, 
which is not the case for elements of H'' . One has however: 

Lemma 11. For any s € (1/2,2], the map r is bounded from the subspace 0/ H" consisting of 
functions that vanish at and n into Hq . 

Proof. Instead of considering the restriction operator r as before, we are going to consider the 
operator f defined on continuous functions as 

(r7)(t) = /(t)-^(/(0)(vr-t) + /(l)t). 

Note that ff = rf if /(O) = /(tt) = 0, so that the statement will be implied by the fact that f is 
shown to be a bounded operator from H'' to Hg. Therefore, instead of considering A as before, 
we consider the operator A — (Aq^^tA"*/^ — r), which has matrix elements 



A — 



-^m^ " m + n odd, 
m + n even. 

^ (m""^?!^"* — 1) m + n odd, 
m + n even. 

Note that s = 2 is a special case since one then has A = as a consequence of the relation 
Aof = rA. 

In this regime, we have as before < C||A||, but this time A is defined via its matrix 
elements by 



Ari 



a m > n. 



Computing {A(pm, Aipn) for m > n as before, we note that there is a difference between the case 
s < 1 and the case s > 1. We obtain: 



\{Alfi„^,Alfin) \ < 



s > 1, 



For s < 1, we now make the choice fm,n — rn'' ^ri'^ ^ , where e > is chosen sufficiently small 
so that 2s — e > 1 (this is always possible since s > ^). With this choice, we obtain 



m>0 



e-2s„2s-l-e 

\ I / 7 lit. If. -r- I / 7 lit. 

fn 



and similarly for the other term. This calculation also works for the case s = 1, so that the case 
s > 1 can be obtained in an identical manner (set for example e — ■ 

Consider now the operator Ca given by (^£af){t) = dff, endowed with the boundary condi- 
tions /(O) = fin) = and /"(O) = -a/'(0), /"(tt) = afin). Since the domain of the square of 
the Dirichlet Laplacian is X>(A^) = {/ G V{Ao) \ Aof € V{Ao)} = V{Cq) we have Co = A^. 
The following lemma shows that Ca for a 7^ can still be viewed as a perturbation of Ag. 
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Proposition 12. Fix a G M and e > be arbitrary and define the linear operator A: Hq 
Ho by 

Af = I'm', - /'(7r)<5;. 

Then, the operator Ca = Ag + a A is the generator of an analytic semigroup on Hq . Furthermore, 
this semigroup coincides with the one generated by Ca, so that Ca = Ca- As a consequence, we 
obtain the identities Ti" = Hq" for every G (— f , f )• 

Proof. Note first that A is well-defined since it follows from standard Sobolev embedding theor- 
ems that /' is continuous for every / e H5+^ and therefore for every / e Hq^^. Thus, S'q and 

S'^ can be considered as elements of the Sobolev space Hq ^ ^ . Since Aq generates an analytic 
semigroup on Hq for any a G K, it follows from applying |Hai09[ Proposition 4.42] once with 
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B = Hg ^ and once with S = Hg ^ that Ca is the generator of an analytic semigroup on Hg for 
every a G (— |, |) and that the corresponding scale of interpolation spaces satisfies V." = Hg" 
for every a € (-§, |). 

It therefore remains to show that the semigroup Sr generated by Ca coincides with the 
semigroup Sr generated by Ca- Since, for any u S Hg, any r > and any t G (0, tt), we have the 
identities 

drSu{T,t) = ~dfSu{T,t), drSu{T,t) = -dtSu{T,t), 

it suffices to show that SrU £ 'D{Ca) for r > 0. Since we already know that = Hj for 
example, we have (S'ru)(0) = {Sru)(n) — 0, so that only the second set of boundary conditions 
needs to be checked. For this, writing for the semigroup generated by Ag, note that we have 
the identity 

SrU = S°u + a[ S^_^ASrudr + a [ S'^_^A{SrU ^ Sru) dr- (12) 
Jo Jo 

Therefore, the first term in this equation belongs to Hg. Furthermore, it follows from the 
definition of A that the Hq norm of the third term is bounded by C Jq {t — r)~^~^~'^\\SrU — 
Sru\\ii2 dr- Since we know already that Hg ='Hi , and since St-u S for every a > 0, it follows 
from standard analytic semigroup theory that jjS'rW — S'ru||jj2 < C|r — t|. So that the third term 
in (fT2)) also belongs to Hg. The second term can be rewritten as 

/ S^_^ASrudr ^ -Aa^ASrU- S° / S°ASrudr- 
Jo Jq 

Collecting all of this, we conclude that we can write 

SrU = —aA'f^'^ASrU + RrU, 

where RrU G Hq. On the other hand, using an approximation argument, one can check that if 
/ e C^, then g = A^^Af satisfies the boundary conditions ^"(0) = /'(O) and g"{'n) = —f'{n), 
from which the claim follows at once. ■ 

Corollary 13. For every a G we have the identity Ha — H^" ([0, tt], R'') . For every 

a e (i,i], we have the identity Ha ^ H^" ([0, tt], R'') n Co([0,7r],R''), where Co([0,7r],R'^) 
denotes the set of continuous functions vanishing at their endpoints. For every a G 
we have Tia = H^" ([0, vr], R'')/ ^, where the relation ^ identifies distributions that differ only 
by a linear combination of Sq and 6,^ - 

Proof- By Proposition [T^l we already know that Ha = Hq" for a € [0, ^]. The claim for a < ^ 
then follows from Lemma [10] while the claim for a G (1; ^) follows from Lemma [11] The 
remaining claims follow from duality. ■ 



3.3 Well-Behaved Projection Operators 

We will later identify the stationary distribution of the SPDE ((H]) by using a finite dimensional 
approximation argument. When projecting the equation to a finite dimensional subspace, the 
most natural choice of a projection would be to use the orthogonal projection H„ onto the space 
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spanned by the first n eigenfunctions of C, but it transpires that these projections don't possess 
enough regularity. Instead, we wih need to use, in some places, the operators n„ given by 

ilnX = y^^- — -{x,ek)ek, (13) 
fc=i 

where the are the eigenfunctions of C. The purpose of this section is to prove the required 
regularity properties for n„. 
We use the Holder norms 

\x{t) — x{s)\ 
» + ||ii|oo + sup — — if a; e and 

|a;||ci+» = S 

else 

where a G [0, 1) and write 0^+" = {x e \ \\x\\ci+^ < oo} and = {x e 0^+" | a;(0) = 

x{Tr) ^ 0}. 

Lemma 14. Let fk : [0, tt] ^ R 6e defined by fk{t) = sin(fci). Define n°a; = LLi ^(^^ fk)fk- 

a) Let Fn be the Fejer kernel given by Fn{t) = i(sin(^)/ sin(|)) for all t e [— 7r,7r]. Then 
tl'^x — ~jFn * X, where x is the antisymmetric continuation of x. 

b) ||n°j:||ci+» < 27r||x||ci+c for all x G 0^+" and all a G (0, 1). 

Proof, a) Since x[t) sin(fct) dt — 2{x, fk) and J^^ x{t) cos(fct) dt = 0, it follows from trigono- 
metric identities that ^ 

{x, fk)fk{s) cos(fc(s - t)) dt. 

The result then follows from the fact that Fn{t) — 2X]fc=i ^'^^ cos kt. 

b) This follows directly from part a, using the definition of the C^^"-norm and properties of 
the convolution operator. ■ 

Lemma 15. Let a G (0, 1) and let x G G^^"' with x(0) — x{Tr) — 0. Then there exists a constant 
c > such that the bounds 

1) ||Mlci+= <cfci+", 

2) ||efc-/fe||ci+» <cfc", 

3) \{xjk)\<c\\x\\ci^.k-^"^, and 

4) \{x,ek-fk)\ < c||x||ci+=fc-2— 
hold for every /c G N. 

Proof. The first bound is standard. The second bound follows immediately from lemmalSJ part b. 
For the third bound, we use partial integration to get 

{x,fk) = T x(t) cos( kt) dt = -y^{-iy ±(t+^Tr)cos(ks)ds. 
k Jo Jo k ' 

Writing = sup5_^( ■^^^^^r^^s^; it is easy to see that each term of the sum is of order 0(y\x\a/k^'^"^ 
and the claim follows from this. 

The bound on {x, e^ — fk) follows similarly: if g is any C^^" function with 17(0) = (?(7r) = 0, 
we can use integration by parts to get 

1 

g{t) sin(A:t) dt — — I g{t) cos(kt) dt, 
k Jo 

g{t)e-'''dt^^J^ me-^'dt, 

and similar results for integrals against cos kt and e^'"'^'^"*' . As above, these expressions are 
bounded by 0{k~^~°'). The claim now follows from lemma |9l part b, by absorbing the slowly 
(7) • 

varying terms gj: into g. ■ 
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The following lemma collects all the properties we will require for the operators li„. These 
will be used in the proof of propositions I26[ below. 

Lemma 16. Let (e„)„gN be an orthonormal system of eigenf unctions of C and denote by n„ 
the orthogonal projection of % onto En = spanjei, . . . , e„}. Define n„ : % — > £"„ as in (jl3p . 
Then the following statements hold: 

a) n„ o n„ = n„. 

h) n„a; ~> X in Ha as n ^ oo for all a g M. 

c) ||n,i||w„ < 1 for all n eN and a €R. 

d) For every 0<a</3<l we have ||n„||(^i+;9_^(^i+Q < cxi. 

e) Let Q < a < P < 1/2 and x e . Then ||n„x — — >■ as n oo. 
Proof. Statement a is clear from the definition of n„. Let x — X]fc°=i ^k^k G Ha- Then 

n ^ oo 
X - tlnX = ^ -XkCk + ^ XkCk 
k=l k=n+l 

and thus, writing for the eigenvalues of —C, 



|a; — Il„a;||^^ — ||(— £)"(a; — n„x)||j^2 — ^ — A^"x^ + ^ Xf." 



oo 



fe=l fc=n+l 



as n — oo. This proves statement b. Similarly, we have 

l|nna;||^^ = ^ — ^ Xl"xl < ^ A^"x| = 

fc=l fc=l 

which is statement c. 

d) From lemma [TS] we get |lefc||ci+Q < ||/fc||ci+° + ll^fc — /fc||ci+° ^ ck^'^". Using this and 
the other bounds from lemma [T5l we obtain 

" n — k 

||n°a;- na:||ci+= = ^ \^{x,ek - fk)ek + {x,fk){ek - /fc)||ci+„ 

k=l 



oo 

- - /fe)|||efe||ci+= + |(a;,/fe)||lefe - /fc||ci+°) 

k=l 

oo 

<C||x||ci+.^fc-^- 



A,-l 

oo 

^,-l-(/3-a 

Since we already know that 11^ satisfies the requested bound, the claim follows. 

e) Let £ > 0. We can write x G Cj'*''' a,s x — y + z with ||2/||ci+« < £ and z e with 
z(0) = z(T) = 0. This gives 

1 1 n-n*^ 1 1 C^+" ~ ^ 1 1 *^ *^ 1 1 C^+" ' 

Because 1 + a + 1/2 < 2, we have ||z||ci+q < c||z||h2 for all z G H^. Corollary [T3] gives 
n Co = and thus 

M II O ^ II II rti/2 

as n — ?> CXI by part b. Since we can choose e > arbitrarily small, this completes the proof. ■ 

4 The Linear Case 

This section gives the proof of theorem H] for the linear case / = 0. 

Lemma 17. Let C and x be given by definition\^ r/ien Qq'^"'"^'^"*" = A/'(a;, — . 
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Proof. Since, for / = 0, the components of the solution of ([H]) are independent, it suffices to 
work in dimension d — 1. First consider the unconditioned process, described by ([5]). It is easy 
to check that p satisfies 



p{t) = e-*/2™p(0) + \j\f^ Q-ii-r)/^^n 

and thus 

g(t)=x_+2(l-e-*/2™)p(0) + y^ / [ e-(''-'')/2" dw(r) dt;. 
The mean of this process is 



JO 



Xo{t)=E{q{t)) =x_ 
and, since p(0) is independent of w, the covariance function can be found as 

Co(5,t) = Cov(g(s),g(t)) 



m I rs rt ruAv 



^0 JO 



for all s, i € [0, T]. Evaluating the integrals and combining the resulting terms allows to simplify 
this to 

Co(s, t) = 2(s A <) + 2m{e~''^'^ + e-*/^" _ ^-\s~t\/2m _ 

Denote the mean and covariance function of the process conditioned on q{T) = xj^ by x 
and C, respectively. From |HSVW05l equations (3.15) and (3.16)] we know that 

x{t) = xo(t) + Co(t, T)Co(T, T)-i(a;+ - x^) 

and 

C{s,t) = Co(.s,0 - Co(s,T)Co(T,T)-iCo(T,t) 
for all s,t £ [0,T]. The covariance operator of (5g'^"'"'"''^+ is then given by 

c/(5)= r c{s,t)f{t)dt. 



To complete the proof we have to verify that x and C have the required form. The following 
facts are easily checked: 

(i) the first the derivatives Co{s,t), 9(Co(s,i) and dfCo{s,t) are continuous at t = s and the 
third derivative at t — s jumps according to 

(ii) the derivative boundary conditions 

2md^Cois,0) = dtCois,0), 2md^Cois,T) = ~dtCo{s,T) 

are satisfied; 

(iii) the left boundary condition 

holds; and 

(iv) LCo{T,t)^0. 

Clearly, by (ii) and (iii), the mean x satisfies the required boundary conditions ([5]) and by (iv) 
it also satisfies Lx — 0. From the definition of L and in particular from properties (i) and (iv) 
we can deduce 

-LC{s,t) ^S{t~s) 

and using (ii) and (in) we deduce that C{t,s) satisfies the boundary conditions ([7]). Thus C is 
the Green's function of —C and we can deduce that C = —C~^ as required. ■ 
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Proposition 18. Consider the L"^ (^[0,T],M.'^'j -valued equation 



dy{T) = C{y{T) - x) dr + v^dwir), y(0) = yo (14) 

where yo €E L'^([0, T], M'') . Then the following statements hold: 

a) Equation (jl4p has a unique, global, continuous mild solution. 

b) For every a < 3/8 the solution y is a.s. continuous with values in Ha- 

c) The distribution Qg'^"' is the unique stationary distribution for 

Proof. From lemma [T7I we know Q^'^^'^'^^ = Af[x, {—C}^^) . Thus we can apply [HSVWOSl 
lemma 2.2] to get that (|14l) has a continuous, L^([0,T],M'^)-valued mild solution and ly is its 
unique stationary distribution. 

Let Xk, fc e N be the eigenvalues of — £. Then, using lemmaHl tr(— £)~^'^ = SfceN ^k'^^ < °° 
if and only if /? > 1/8. Thus, for example by applying |Hai09, Theorem 5.13] to y — x, the solution 
takes values in Ha for every a < 1/2 — 1/8 = 3/8 and is continuous by |Hai091 Theorem 5.17] 
(see also |DPZ92) for very similar results). This completes the proof. ■ 

The regularity of the solution given in proposition 1181 is consistent with the regularity of the 
target distribution Qq'^"'"^'^^: the process i in ([5]) is continuous and lives in H^/^~'^ and thus x 
is in H^/^""^ for all e > 0. On the other hand, corollary [T3] shows that 'Ha C H^" and thus that y 
also takes values in H^/^~^ for all £ > 0. The following lemma provides an additional regularity 
result for x in stationarity. 

Lemma 19. Let a < 1/2. Then x e C^+" for Qq'^''^''^^ -almost all x. 

Proof. The result is a direct consequence of |Hai09[ Corollary 3.22]: Let ek be the eigenfunctions 
of —C with corresponding eigenvalues A^. By lemma [T71 part b, the random variable 

oo 



fc=l ^ 

where the 77^ are i.i.d. standard Gaussian random variables, has distribution Qg' ^' . We 
have to show that the derivative 



00 



X' + x' = Y. 



Vk 



is a-Holder continuous. 

Let 5 e (2a, 1). By lemma[9]we have lle^^lU = 0[k), \\e'l\\oo = 0{k^) and Afc = ck^ + 0{k^) 
for some c > 0. This gives 

oo / oo / / oo 

o2 V^ll '^fc l|2 r,2 V^ll et, ||2-5^. / Sk \^ ^ Sr^ ^ 



fe=l ^ " fc=l ^ ^ k=l 

Thus the conditions of jHai091 Corollary 3.22] are satisfied, and we get the required Holder 
continuity. ■ 



5 The Nonlinear Case 

In this section we complete the proof of theorem |4l The proof is split in a sequence of results 
which identify the target distribution Q j' " ' , determine the regularity properties of the 
drift Af, give existence of global solutions to the SDE and finally identify the stationary 
distribution of this equation. 

Lemma 20. Assume that /: M'' — !■ R'^ is such that the SDE ^ a.s. has a solution up to time T. 
Let II = Qj^^' and v = Qg'^"' be the distributions on L^([0,T],R'') from definitions^ 
Then the density ip — ^ is given by 

fi^) = exp(m{f{x+), x{T)) - m{f{x^), i(0)) 
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^ m{Df{x{t))x{t),x{t)} - {f{x{t)),x{t)) + llfi^m' dt) 
where Df is the Jacobian of f and Z is the required normalisation constant. 

~ Ox . . ... . . ^ Ox 

Proof. Let fix = Pj' ~ be the unconditioned distribution of i in ([6| and let S'x = Pq ^ the same 
distribution, but for / = 0. Then the Girsanov formula, e.g. in the form of [HSV07[ lemma 9], 
gives the density of fix w.r.t. t'j,: 

dfLx I . 



-^(i)=exp(/ {f{x{t)),dx(t)) + :, {f{x{t)),-xit)-f{xit)))dt) 
av± \Jq ^ jo III' ' 

where a; is a deterministic function of x via the relation xif) = + i(s) ds. Since 1 1— > f(x(tf)^ 
has bounded variation, we can use partial integration to get 

(/(x(t)) , dx(t)) = ifi^iT)) , i{T)) - {f{x{0)) , m) - r {x{t),Df{x{t))x{t)) dt. 

Jo 

Substituting this expression into the formula for dfL±/dD± and using substitution to switch from 
X to X gives 

^(x) =exp(m(/(x(r)),i(T)) -m(/(:r(0)),i(0)) 

-m^ {x{t),Df{x{t))x{t))dt+\j^ {f{xit)),xit))^\fixm'dt) 

~ Ox. . . ... . ^ Ox. 

where fix — Q f ^ is the unconditioned distribution of a; in (|6]) and Vx = Qq ^ is the correspond- 
ing distribution for f — 0. Now we can condition on x{T) — x+, e.g. using [HSV07I lemma 5.3] 
to get the result. ■ 



Lemma 21. Let a € [0, 1). Then there is a c> such that 

lci+° 



|i|lL+"<c||x||^||x|| 



for all X e Ci([0,r],] 

Proof. The claim for a = is trivial, so we can assume a 7^ 0. Assume first the case d = I. 

\x{t) — x[ 



Write \x\a = s\xps^t ^^'"\t-s\^J'^^ let t e [0,r] such that \x{t)\ = Then 



^ \x{t) - x{s)\ ^ \\x\\oo ~ \x{s)\ 
' \t-s\^ - \t-s\^ 

and thus |i(s)| > ||i||oo — for all s e [0,r]. This allows to conclude that > 
on an interval of length at least TA ||i;||i!,^"/(2|i;|Q)"'^/". Since we assumed d=l, this gives 

11 llill^" \ . /T 



2i/-|i;|y"/ V4" '22+i/-|i|^/^ 

and by solving this inequality for ||a;||oo we find 

<c||x||^max(|i:U,||x||oo) <c||x||^||a;||ci+. 

for some constant c. 

For d > 1 we apply the inequality componentwise: Since, for z e R"*, we have ||z||2/\/rf < 
\\z\\oo < \\z\\2, we get 

||i|lL+" < c max (||x,||^(||x,||oo + ||ij||oo + l^:,^)) 

j = l,...M\ / 

<c||a;||^(||x||oo + ||i||oo + |i;|a) 
where c is increased as needed. This completes the proof. ■ 
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The following bound for the density will be used to show that the stationary distributions 
of approximations for the sampling SPDE Q are uniformly integrable. 

Lemma 22. Let Lp he the density from lemma \2(K U ~ log 93 and v — Qq'^"'^'^+ and a £ (0, 1). 
Then for every e > there is an M > such that for v-almost all x we have 

U{x) < e\\x\\li+^+M. 

Proof. We bound the five terms in U one by one. For simplicity we denote all constants in 
the following estimates by the symbol c, the meaning of which changes from expression to 
expression. 

1) Using the Cauchy-Schwarz inequality we get the bound 

{f{x+),x{T))<\fix+)\\\x\\^<\f{x+)\\\x\\ci+^ < e\\x\\l^+^ + c. 



2) A very similar argument gives —{f{x-),x{0)) < £||a;||?,i 



c. 



3) We can use Young's inequality together with lemma [H] to conclude that for every e > 
there is a c > such that 

||i||oo < £|l2:||ci+° + c||a;||oo 
for aU x e L2([0,T],M'^). Thus we have 

{Df{x{t))xit),x{t))dt<\\j^f{x)\\^ \\x\\^ 
<T\\^f{x)\\ \\x\\ 

- II ^ ^lloo II lloo 



< 



(*^ll/(a;)||ci+° +c||/(a;)||oo) (£||a:^||ci+= +c||a;||oo)- 



Since / is differentiable with bounded derivatives we have ||/(a::)||ci+a < c||x||ci+a + c and 
by assumption there is a /? < 1 such that < + c. Using these estimates we find 



T 

(D 







f{x{t))x{t), x{t)) dt < [e\\x\\c^+. + c\\x\\i + c) {e\\x\\c^+. + c\\x\\ 



and thus for every £ > there is a c > such that 







{Df{x{t))x{t),x{t)) dt < e\\x\\l,+. + eMl + c\\x\\l^P + a 



Since /3 < 1, this gives the required bound. 

4) Using the Cauchy-Schwarz inequality again, we get in a similar way 

T 

{f{x{t)),x{t)) dt < ||/(a;)||2||i||2 < c(||x||^ + c)\\x\\^ < eMl^+^ + c 



5) Finally, we have - \f{x{t))\'^ dt < 0. 

Combining these bounds gives the required result. ■ 

Lemma 23. The drift Af defined by pop is locally Lipschitz from TLi/n to H-r/ie- Furthermore, 
one can write M = Mi + M2 + M3 such that A/i does not depend on x and such that the bounds 

mx)\\n_y,,<c{i + \\x\\i^^j, m) 

\\M2ix)-M2{y)\\n^,,,, <c\\x~ y\\n,,AMn,,, + \\y\\n,„), W 
||A/-3(a:)-A/-3(2/)||«_3/,, <c||x-2;||«,,,(||x||„,,, + ||2/||„,^j' (USt) 

hold for all pairs x,y G TLi/i and for some constant c > 0. 
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Proof. We use the characterisation of the spaces Ha from coroUarv [T51 and in particular the fact 
that if x G 'Hi/i, then x also belongs to H^. Since we assumed that /_,- and its derivatives up to 
the second order are globally Lipschitz, this implies that fj{x), difj{x) and dijfk{x) all belong 
to and their norms are bounded by multiples of that of x. 
Let now x G H^. Then the following statements hold: 

• fi{x)dkfi{x) G since is stable under composition with smooth functions, 

• dtXi(^difk{x) — dkfi{x)) G L^, for the same reason 

• dtXi dtXj dfjfk{x) G H* for aU < < — 1/2 since in this case C H* by Sobolev embedding, 

• dfxj(^djfk{x) + dkfj{x)) G H^^ since H^^ is stable under multiplication by H^-functions, 

• fk{x^)dtSa G H* and fkix+)dt6T G H* for every t < -3/2. 

It follows that JV maps H1/4 into Ha for every a < — |. In particular, it maps H1/4 into 
as stated and the bound p5b|) holds. We then define A/i as the term proportional to 
fk{x-)dtSo — fk{x+)dt5T, As as the term proportional to dtXi dtXj dfjfk{x) and Af2 as the sum 
of the remaining terms in the nonlinearity. With these definitions at hand, the bounds (jl5bp 
and (jl5bp follow easily. ■ 



Proposition 24. For every initial condition xo G L^([0, T'],R'^), the stochastic evolution equa- 
tion (jn]) has a unique maximal local solution (x, t*). The solution satisfies x{t) G H1/4 for every 
T < T* a.s. and sup^^^, ||a;(T)||L2 = 00 a.s. on the set {t* < 00}. 

Proof. Define 

g{T) ^ S{T)xn + V2 I S{t - cr) dw{cr). 
Jo 

Let R> U > 0. For x: (0, U] Hi/i continuous, define 

11x11.= sup rV4||^(^)|| 

re(0,(7] 

and let X be the space of all such x with a;(0) — g{0) and ||a;||, < 00. Then (^X, || • ||») is a 
Banach space. We find 



||fflU< sup T'/^(^\\xoh,+V2\\ r SiT-a)dw{a)\\ ) 
< ||xo||l2 + \/2i? sup S{t - a)dw{a)\\^ 



=: IIxoIIl^ + Cr 

and thus g E X for every U < R. Define a map Mg : X X hy 

Mgx{T)^ f S{T-a)Af{x„)dcr + g{T) VrG[0,C/]. 



By the definition of a mild solution, local solutions up to time U coincide with the fixed points 
of this map. 

Let B{g, 1) C A" denote the closed ball around g with radius 1. By lemma the nonlinearity 
A/": Hi/4 — > H-T/iQ is locally Lipschitz and thus for all x,y E B{g, 1) we have 



Mgx-Mgy\\^< sup t1/4 / \\S{t - <j){Af{x,) - M{ya))\L da 

rG(0,(7] -'1 ^'^ 



M2{Xa) - U2{y^)\L \W^{Xa) ~ Mz{ya 



(t)9/16 (t-(t)7/1 



da 



' ^ I II "^cr \ \u ( ' 

Jo " ll«l/4V 



(r-cr)9/16 + (T-cr)7/16 



2 



< sup cr / / - 2/^11^, ( -^^^7^^^ + ] da 

re{o,u] 
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<cU'/''\\x-ymi + \\x\\, + \\y\Uf, 

where c changes from hne to hne. Similarly, we have 

||AV - g\U < cU'/'' Ml < cU^I^^ (^\x ~ g|U + |l.g|U)'. 

By choosing the final time U sufficiently small, we can then make sure that Mg is a contraction 
on the ball -B(g, 1) and, by the Banach fixed point theorem, Mg has a unique fixed point. This 
gives a unique local solution of © up to time U . 

By iterating this procedure, every time starting with the final point of the previously con- 
structed segment, we obtain a solution up to a maximal time r* < R. Since the length of each 
segment of this solution only depends on the L^-norm of its starting point, we see that t* < R 
implies sup^^^. ||a;(T)||L2 = oo. Taking R oo completes the proof. ■ 

Even if / is globally Lipschitz, the dtXi dtxj df^fk-teim causes the nonlinearity J\f to be only 
locally Lipschitz. Thus, showing the existence of global solutions to the SDE ^ will need some 
care. 

Proposition 25. For every initial condition Xq G L^([0,T],R'') the SPDE ^ has a unique 
global solution. For every r > the solution satisfies E(||2;(r)||^2) < oo. 

Proof. From proposition[M]we know that ^ has a local solution (x, Tmax)- Let y be the solution 
of the linear SPDE from proposition [TSl i.e. 

y{T) ~ S{t) (xq — x) + a/2 / S{t — a) dw{a) + x, 

Jo 

and define z{t) — x{t) — ^(t) for every r e [0,Tiiiax)- Then z satisfies the stochastic evolution 
equation 

dz{T) = Cz{t) dT + N{z{t) + y{T)) dT, z(0) = 0. 
Thus ||^;(t)||^2 satisfies 

= 2{z{tICz{t) +N{z{t) + y{T))) 

= -Am^{^lz{T),^^^z{T)) ^ Am\dtz{Q)\^ - Am\dtz{l)\^ (16) 

- {dtz{T),dtz{T)) + 2{z{t)M{z{t) + y{T))) 
<-c\\z{T)\\l,+2{z{T),N{z{r)+y{T))) 

for some c > 0. This formal calculation can be made rigorous by a standard approximation 
argument, using for example Galerkin approximations. 

We require a priori bounds of the form (^z,N{z + y)) < c\\z\\^2 + s\\z\\^2 + c where e > is 
small enough to be compensated by the negative ||z|jy2-term in ([TC)) . 

In order to obtain the required bounds we consider the five terms from the definition of Af 
individually. For the purpose of these estimates we denote all numerical constants by c > and 
only track the y-dependency of the bounds explicitly. For the first term we get 

(zk, -Mz + y) dkMz + y)) < c||z||l2||/(z + y)h2 < c\\z\\l2 + cMl^ + c. 

For the second term we find 

(zk, dt(z, + yi) dt(zj + yj) dfjfk{z + y)) 

ZkdtZidt{d^fkiz + y))dt+ / Zk dtyt dt^Zj + yj) dfjfkiz + y) dt 



1 /.I 

,2 



dtZk dtZi d^fk {z + y)dt- / Zk z^ difk {z + y) dt 

Jo 

+ Zkdtyidt{zj + yj ) df^ fk{z + y)dt 
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< c||z||hi + c||z||l2||z||h2 + c||z||l°o Hyllni (II^IIhi + WvWr^)- 
For the third term we have 

{zk,-dt{z, + y,) (5Jfc - dkfi)) < c||z||l2 (||z||hi + ||y||Hi)- 
The fourth term can be bounded as 

{zk,^Uz^+y^){^^.fk+^kf^)) 

= {zk,dfzi {djk + dkfi)) 

»i .1 

dtZkdtVr [difk + dkfi) dt- ZkdtVi dt{zj + y^) {df^fk + dfkfi) dt 

Jo 

< c||2:||l2||z||h2 + c||z||Hi||y||Hi + c||2:||L-||y||Hi(lk||Hi + l|y||Hi)- 

Finally, for the fifth term, involving the derivatives of Dirac distributions, we get 

{z,dt6o) = -z'(0) < c||z||h^, {z,-dtSi) = z'{l) < c||z||h" 
for every a > 3/2. 

To convert the bounds into the required form first note that for every s G (0, 2) the interpol- 
ation inequality (see e.g. [Hai09[ corollary 6.11]) gives \\z\\^s < \\z\\'^'^\\z\\^2 and, using Young's 
inequality, we can for every e > find a c > such that 

Using this relation we find a c > such that 

IkllL^llz/ll^i < + \\\y\\j,. < c\\z\\l. + cMti^ < 44h + e\\z\\l. + cMt,.. 

The terms of the form ||z||hi llyllni can be bounded using the relation 

||2||l-||z||hi < c||z||h3/4||z||hi < c ||z||l2^||z||^''/ ||z||l2^||2;||h2^ = c\\z\\%^\\z\\K2. 
Applying Young's inequality with p — 16/7 and q — 16/9 we find a c > such that 

ll^llL-lkllHilbllHi < c||z||^/* ||z||L//||y||Hi < elklln^ +c||z||l2||2/|Ih/^ 
Combining all these estimates, we find that for every £ > there is a c > such that 
(z(r), AA(z(r) + y{r))) < c(l + \\y\\]^T) \\z\\l + s\\z\\l,. + c(l + ||y||^, + ||y||4 ^) 
and substituting this bound into ([T5| for small enough e > we get 

< c(l + \\y\\T)Hr)\\l + c(l + 11,11^2 + Ml,.). 

Gronwall's inequality gives 

Mr)\\h < c r(l + \\y{a)\\]^/') Hi + \\yir)\\l + \\yir)\\l,r) dr 



exp 



{l + \\y{r)\\'^l')dr)dcr (17) 



[i + \\v{^)\\i-^ + \\y{<y)\\h)d<y. 

Thus II 2 II cannot explode in finite time and from proposition [53] we get Tmax = oo. 

By proposition [TH] we have y G ([0, t], H^) . Hence, by Fernique's theorem (see e.g. [HaiOQi 
theorem 3.11]), 

exp(e / ||y(f )||hi ^'") 1 < oo 



for sufficiently small e > 0. Thus, using the fact that 16/9 < 2, we see that the right hand side 
of (fT7|) has finite expectation for all t > 0. ■ 
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Now the only part of theorem 2] which we stih need to prove is the statement about the 
stationary distribution of This can be done using a finite dimensional approximation argu- 
ment, similar to the proofs in [Zab88' and jHSVOTI section 3]. Since these articles assumed that 
U was bounded from above and also assumed different regularity properties for the drift, the 
proof needs to be adapted for the situation here; to allow for easier reading, we include the full 
argument instead of just enumerating the required changes. 



Proposition 26. The distribution Q 



invariant for ^ . 



X 'T x-i- X 'T x-i- • 

Proof. Let (p be the density of /i = Qj' ' + w.r.t. u = ^' as given by lemma [201 and 
let U — logiy9. Then we can compute the derivative oi U at x £ T-ii/i in direction h G "Hy/ig as 

{DU{x),h) = 7n fk{x+)hk{T) ~ m fk{x^)hk{0) 

+ 2 j (-fidkfi + mxiXjdfjfk - ^Xi{difk - dkfi) +mxi{difk + dkfS)hk{t) dt 

^ {N{x),h). 

Here we used the fact that, by corollarvfT51 h € 'Hy/ig implies /i(0) = h{T) = 0. This shows that 
the function U is Frechet-differentiable with derivative M . Let n„ and II„ be as in lemma [16] 
and define the approximations 

AA„ = (t/on„)' = n„AA(n„-) 



for n e N. 

Consider the n-dimensional SDEs 



V2n„ dw(r), y„(0) = linXo 



and 



dxnir) = £x„(r) dr +7V„(a;„(r)) dr + \/2n„ (iw(T), 2;„(0) = n„a;o. 



Then, by finite dimensional results, the stationary distributions i^n and Hn of 2/„ and x„ respect- 
ively are given by 

-1 _„j ^Mn 



o n„ 



and 



dvn 



exp(t/on„). 



Define the semigroup (7';')r>o on Cb('H,M) by Vl'ipix) = E^{ip{xniT))) for all x € and 
ip G Cb(H,R). Since the process Xn is ^„-reversible, we have 



■H 



(p{x)r^ij{x)dHn{x) ^ / tp{x)v:;(p{x)dfin{x). 



(18) 



n 



for every ip,iP € Cb(n,R). 

We need to find the limit of (jl8p as n — > oo. For this, we first show that Xn — ?> a; in H1/4, 
uniformly on bounded time intervals: Let U > 0. Then we have 



\xn{T) - x{t)\\^^^^ < II (n„ - /) (^(r)xo + V2 S{T-a) diy(a)) 



S{t - cr){j\fnix{a)) - Af{x{a))) da 



Sir - a){Nn{xn{(T)) - Nn{x{a))) da 



«l/4 



for ah T G [0,[/]. 

From the definition of || • Wn^ and the asymptotics of the eigenvalues of C in lemma[S]we get, 
for any /? > a, that there is a c > such that the bound 
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holds for all X E Hp and all n e N. Let /3 e (1/4, 3/8). Then we know from proposition [T51 that 
T H> S{t)xo + \/2 S{t — a) dW{a) is a continuous map from [0, U] into Hp. Combining these 
two statements, we find supQ<^<[j hiT) — >■ as n — > oo. 

From lemma [23] we know that Af is locally Lipschitz from H1/4 to H^t/iq. By lemma fTHl 
part c, there is then a constant Kr > such that 

\\j^n{x) ~ Afn{y)\\^_^^^^ < Kr\\x-y\\n,^, 

for all n e N and all a; and y with ||?/||-Hi/4 < f- Thus, the Afn are also locally Lipschitz. 

We can find p,q > I such that p ■ ^ < 1 and 1/p +l/q — 1. For I2 we get then 

/2(t) < r\\S{T - a)(A/-„(x(a)) - AA(a;(a))) da 







I ^1/4 



-^d.) [j^ \\U^{x{a))-U{x{a))\\l^_^^^^da) . 

The right hand side is independent of r and converges to as n — >■ 00 by dominated convergence, 
using lemma [TBI part b. 
For n £ N define 

T„^r = infjr € [0, U] | |la:(T)|| > r or |1x„(t)|1 > r}, 
with the convention inf $ = U. For r < Tn^r we have 

h{T)<Krj ||5(T-cr)|| ||x„(cr) - a;(cr)|| rfcr 

II II H_7/16— >tti/4 II 'llttl/4 



and consequently 



||a;„(T) - a;(T)|| < sup (li{a)+l2{a)) 

Using Gronwall's lemma we can conclude 

||a;„(r)-x(r)|| < sup {li{a) + hia)) ■ exp(cKr [ I, da) 

for all T < r„_r- As we have already seen, the right-hand side converges to as rt — >■ 00. 

Now choose r > big enough such that supo<T-<t/ ||a;(T)|j < r/4. Then for sufficiently large n 
and all T < T„^r we have ||a;„(r) — a;(T)|| < r/4 and thus supQ^^^j^^ ^ ||a;„(T)|| < r/2. This implies 
Tn,r = U for sufficiently large n. Thus we have x„ a; in C([0, U],Hi/i) a.s. 

Let < a < /3 < 1/2. Define the semigroup {Vr)r>o on Cb('H,R) by VrV{x) = ¥.^{ip{x{T))) 
for aU a; e "H and ip G Cb('H, M). Then, by dominated convergence, we have V^ipi^lnx) — Vrfix) 
as n — >■ 00. By lemma [TOl x € C^~^^ for iz-almost all a;. Furthermore, U : C^"*"" — R is continuous 
and thus U{Il„x) — > C/(a;) as n — >■ 00 for (/-almost all x by lemma [T6l part e. 

Finally let c = ||n„||(^i+/3_j_(^i+Q. Using Fernique's theorem we can choose e > such that 

the function exp(e c||a;||^i+^) is z/-integrable. By lemma we can find an M > such that 
U{Ilnx) < e||n„a;||^i+„ +M < ec\\x\\'^i+fi +M for all n e N and i^-almost all x. Then dominated 
convergence gives 

lim / if{x)r:^ipix)dnnix) = lim / ^(n„a;) p;'V(n„x) e^(n"^)dz/(a;) 

ip{x)VrHx)Q^^''^dv{x) = [ ip{x)Vri^{x)dlJi{x) 
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and using ([T^ we get 



n Jn 



Thus the process x is /i- reversible which is the required resuh. ■ 

Propositions [Ml [21] and [211 together imply all claims of theorem [H and so the proof of the 
result is complete. 
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